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■ ABSTRACT 

a 

Aims. In recent years cosmic shear, the weak gravitational lensing effect by the large-scale structure of the Universe, has proven to 
be one of the observational pillars on which the cosmological concordance model is founded. Several cosmic shear statistics have 
. been developed in order to analyze data from surveys. For the covariances of the prevalent second-order measures we present simple 

and handy formulae, valid under the assumptions of Gaussian density fluctuations and a simple survey geometry. We also formulate 
these results in the context of shear tomography, i.e. the inclusion of redshift information, and generalize them to arbitrary data field 
geometries. 

■ Methods. We define estimators for the E- and B-mode projected power spectra and show them to be unbiased in the case of Gaussianity 
Oh' and a simple survey geometry. From the covariance of these estimators we demonstrate how to derive covariances of arbitrary com- 
binations of second-order cosmic shear measures. We then recalculate the power spectrum covariance for general survey geometries 
and examine the bias thereby introduced on the estimators for exemplary configurations. 

Results. Our results for the covariances are considerably simpler than and analytically shown to be equivalent to the real-space ap- 
proach presented in the first paper of this series. We find good agreement with other numerical evaluations and confirm the general 
properties of the covariance matrices. The studies of the specific survey configurations suggest that our simplified covariances may be 
employed for realistic survey geometries to good approximation. 
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1. Introduction 



Despite the fast progress in cosmology during the past years (e.g. ISpergel et al.ll2007 ', and references therein), still only little is 
OO known about those 95% of the energy content of the Universe which is non-baryonic. Especially for the dominant dark energy, 
astronomical observations provide the only source of information to reveal more about its properties. One of the most promising 
methods to shed light on the dark ingredients of the Universe is cosmic shear, the gravitational lensing of distant galaxies by the 
large-scale matter structure. 

Observations of a cosmic shear signal are challenging since the shear, i.e. the distortion of galaxy images by gravitational lensing, 
is small compared to the intrinsic galaxy ellipticities, of order 1 %. Besides, measurements of galaxy ellipticities are further hampered 
by atmospheric seeing, telescope optics and pixelizati on of the im age. Therefore, although theoretical work on light propagation 
in an inhomogeneous universe exists for f our deca des ('Gunn"1967), the first detection of cos mic shear was only reported in 2000, 
independently by Bac on et all (|2000), Kai ser et all (i2000), van Waerbeke et al. (2000) and WittmanetaD (l2000l) . Recently, cosmic 
shear o bservations have re ached statistical errors on cosmologic a l parameters com patible to t hose of other e stablished methods, 
see e.g. Jarvis et"an (l2006h with CTIO da ta, iHoekstra et all (l2()06h . ISemboloni et all (|2006.) and lFu etall (12007 ) with samples from 
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CFHTLS, and lHetterscheidt et al.l (l2007l) with an analysis of the GaBoDS survey. In combination with CMB data, cosmic shear is 



able to break parameter degeneracies because it probes the matter distribution at lower redshifts and on smaller scales than the CMB 
fluctuations. Cosmic shear is also complementary to methods such as type la supernovae and galaxy redshift surveys, its strength 
being that the underlying theory is built on few physical assumptions; in particular there is no dependence on the relation between 
baryons and dark matter. 

The statistical properties of the (Une-of-sight projected) matter distribution measured by cosmic shear are usually characterized 
by the convergence power spectrum P^it), where t denotes the Fourier variable on the sky. If the matter density field is a realization 
of a Gaussian random field, the second-order measure Pk{{) fully describes its properties. This is the case if the density perturbations 
are linear. Equivalently, other second-order statistical measures as the correlation functions, the shear dispersion and the aperture 
mass (for definitions see Sect. |2]i can be used instead of Pk{(). Of all these second-order statistical methods, the shear correlation 
function is most straightforward to measure from data, being insensitive to gaps and the field geometry of the imaged sky region, 
and it provides the full second-order statistical information contained in the data. 

In order to assess strengths and limitations of all these measures it is important to kn ow eiTors and their corre lations on different 
scales characterized by the covariance matrices. Under the assumption of Gaussianity. ISchneider et all (l2002al) . hereafter Paper I, 
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calculated such covariance matrices analytically, starting from an estimator for the correlation functions By further assuming 
a connected survey area, negligible boundary effects due to the finite field size and randomly placed galaxies in the field, the 
ensemble average over galaxy positions of the covariances was taken, and expressions for the covariances as a function of the 
theoretical ^-t were derived. In this work we aim at recalculating the covariance matrices, but now in a Fourier space approach, 
leading to significantly simpler expressions. These are then shown to be equivalent to the result of Paper I. 

In Sect, m we give an over vie w on second-order measures of cosmic shear and their interrelations, using the notation of 
iBartelmann & Schneideil ((2001) and' Schneideij (1200 6'). Moreover, we briefly summarize the E- and B-mode decomposition of the 
shear field. We derive estimators for E- and B-mode power spectra in Sect. [3] and show them to be unbiased under the assumptions 
made in Paper I. The covariance of these estimators and successively the covariances of other second-order measures are calculated 
in Sect. |4] We demonstrate the equivalence of our formulae to those given in Paper I in Sect. |5] In addition, we evaluate the covari- 
ances numerically and discuss our results. In Sect.|6]the calculations of the foregoing sections are generalized to shear tomography, 
i.e. the inclusion of photometric redshift information, as well as to arbitrary survey geometries. We discuss the behavior of our 
estimators for specific examples of survey geometries. Finally, we summarize our findings and conclude in Sect.|7j 



2. Second-order measures of cosmic shear 

In order to obtain information about cosmic shear from a given data set, galaxy ellipticities are measured at positions Gj in the sky. 
The observed complex ellipticity e = ei 4- ie2 contains contributions from the intrinsic source ellipticity and the shear y induced 
by gravitational lensing, i.e. e, = yiOi) + e^. By considering a pair of galaxies at positions & and d^ + O, one defines the tangential and 
cross component of the elUpticity at position ^ for this particular pair as 

ft = -Re(e e"2'^) ; ex = -Im(e e^^i^) , (1) 

where (p is the polar angle of the separation vector (Ukewise for y and e^). Then one defines the shear correlation functions 

UO) := <rtrt>(0) ± <rxrx>(0) . (2) 

Further second-order statistical measures are the shear dispersion defined as the average of the mean shear in circular 

apertures of radius 9, and the dispersion of the aperture mass in circular apertures of radius 6. For an aperture located atO-0, the 
aperture measures are given by 

Map(0) = J Qm rtW ; MAO) = J Qm TxW , (3) 

when {)■ - The weight function Q('&) is axially-symmetric, but otherwise arbitrary. 'Crittenden e t al.l (f 2002b suggested a weight 
function with exponential fall-off. Here Qi§) is chosen to be the same as in . Schneider et al. (2002b, hereafter S02): 

Q(^)^^\^-1^\h(0-&), (4) 
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where H(6 - §-) is the Heaviside step-function. 

The convergence or dimensionless projected surface mass density k is related to the shear via 

^ p (p- (p- 2i0 

y{0) - - dHD(0-&)K{-») with D{0) = ^ ' „ ' . (5) 

In ordinary lens theory a- is a real quantity. If the complex y can be calculated from k, yi and 72 are mutual ly dependent, which l eads 
to a relation that ensures a shear field, generated purely by lensing, to be curl-free as demonstrated in Crittenden et al.l (I2002h and 
S02. As y is a polar quantity, a general shear field can be decomposed in analogy to electromagnetic polarization into a curl-free 
E-mode and a divergence-free B-mode. Although not caused by gravitational lensing, B-modes can be present in cosmic shear data, 
for instance due to the underestimation of noise or uncorrected systematic errors. Further (small) sources of B-modes are source 
redshift clustering, see S02, and the limited validity of the Born approximation, see e.g. lJain et alj (|2000). 

For a pure lensing signal we will set the B-mode contribution to zero later on. Nevertheless we must account for noise in the B- 
mode channel as it is measured by the real-space estimator for the correlation functions used in Paper I, which does not discriminate 
between an E- or B-mode origin of galaxy ellipticities. In order to account for both E- and B-modes, we rewrite the originally purely 
real convergence k{0) as a complex quantity, which reads in Fourier space K{t) - ke{{) + iKsiO- In terms of ke and kb the following 
power spectra are defined via 



(^e(^)%(^')) = (27t)W>(( - n Pe({) , 

{KB(e)r^(e')} = (27:)W\e - 1) Pb({) , (6) 

{kE(e)Kl(n) = (2n)^6'^\e - t) Peb(0 , 

where 5^^\() is the two-dimensional Dirac delta-distribution. Note that all two-dimensional power spectra appearing in this work 
are power spectra of k, which is why the index k is dropped. To lowest order, i.e. in the Born approximation and without lens-lens 
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coupling, the E-mode power spectrum is related to the three-dimensional power spectrum of the density fluctuations via Limber's 
equation 

where is the comoving distance and;\fh the comoving horizon distance. Hence, P{{) is the projection of Pg along the line-of-sight, 
weighted with the lensing efficiency g(x) - dx'p{x')(l -xlx')^ where pix) is the distribution of source galaxies in distance. For 
simplicity we have assumed a spatially flat geometry. From now on we set PEB.{t) = because the cross power is expected to vanish 
for a statistically parity-invariant shear field, as is the case for its corresponding correlation function ^x(ff) - {7i7x}{d)- 
The Fourier transform of (01 reads 

y(0 = e^^'^m , (8) 

where /3 is the polar angle of (, and where k is now a complex quantity. By means of this relation and (|6]l one can derive the 
second-order shear measures as a function of the E- and B-mode power spectra, see e.g. S02, 

Jo 27r Jo 2n 

r d£ t 516Jl{{ff) , r d(( 516Jlm 
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where Jf^ix) is the Bessel function of the first kind of order fi. 

So the above measure s are all linear functiona l of the power spectra, each employing a different filter. As a consequence they are 
interrelated, as shown bv lCrittenden et all ( l2002h . Thus, in practice only one measure has to be estimated from the data, the others 
can then be deduced. As real data contains gaps, masks and CCD defects, the placing of apertures needed for (|yp)(6') and {Mlp}(6) 
is impractical, so that the correlation functions are the preferred primary observables. Due to the broad filter especially of ^+ the 
signal is strong, but mixes information about the power spectrum over a wide range of scales. On the c ontrary, the aperture mass 
statist ics employ a naiTow filter, which can be replaced by a Dirac (5-distribution with an error of only 10% jBartelmann & Schneideij 
Il999h . Thus the signal is small, but probes the power spectrum very locally, and measures on different scales quickly decorrelate 
(see Sect. [Sj. In addition to this, (|9j shows that (M^p)(0) is sensitive only to E-modes, whereas {M\){6), being sensitive only to 
B-modes, can be used to quantify systematic errors etc. The filters of the second-order measures also determine the form of their 
covariance matrices, as will be discussed in Sect. |5] 



3. Power spectrum estimators 

Consider now galaxies at positions 0,- with measured ellipticities e, — eij + ie2,i in a data field covering a solid angle A. Repeating 
a calculation by Kaiser ( 1998), we write 



(10) 



where in the second step e„,, was split up into shear and intrinsic ellipticity. Note that ea{t) is the discrete Fourier transform of the 
real ellipticity component e„ so that e^it) - eai-{) holds. The introduction of n{0) := J]j6^^\6 - 01) regains a continuous Fourier 
transform for the shear. Transforming ya{0) back to Fourier space then yields 

e.(0 = J^r«(nn(^-n + 2</e"'' n{t) ^ ^ d^9 e'Me) ■ (11) 

The explicit dependence of n{0) on the galaxy positions hampers analytical progress. To simplify this expression, we make the 
assumption that the galaxies in the field are uniformly sampled, which means that n{0) can be replaced by its ensemble average over 
all galaxy positions 

E{n{0)) = n H{0) with £ := |~| ^-^ J d^e^j , 

where h - N/A is the mean number density of galaxies in the field and Yl{0) is the aperture function, which is 1 for within the 
field and else. In particular, A - jd^6 Tl(0). With the definition 

A(.):.^. r^e-n(.) (12) 
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one can now write h({) - {2n)^hA((). Note that (2;7r)^A(/') is the Fourier transform of the aperture. If one further assumes that the 
observed field is simply connected and that all relevant angles are considerably smaller than the extent of the field, i.e. |0p A, 
then the aperture function 11 can be approximated by unity, so that consequently A{{) ^ 6^^\t). Through this ( fTTI) takes on the form 

^ nr«W + </ e'^ "' ■ (13) 

In the following calculations terms containing elhpticities to the square will appear. In this case the approximation for Mt) is 
executed as follows 

|A(oi' ^ <J<''(/') m = s'-'-Xe) A(0) = . (14) 

With these considerations at hand we now define the following estimators for the E- and B-mode power spectra, 

^e(^) := ^TTT^ r I cos ip + UO sin2y3|2 - ^ , 
n^AAR{€) Jam 2n 

^TTTa r |-eiWsin2^ + e2(/')cos2/3|2 - ^ . (15) 
n^AAR{() J And) 2n 

These estimators are band powers for disjunct bins I in Fourier space, constructed by averaging over annuli with mean radius { 
and area A«(^) = 2ntM. crl denotes the dispersion of the complex source ellipticities. The Fourier elhpticities e„(^) in ( fT5] l are 
calculated from the data via ( fTOl i. In the following we demonstrate that these estimators are unbiased under the assumptions made 
above. 

On calculating the expectation value, one splits up the image ellipticities into a shear and a source ellipticity part following 
(fTSl l. In order to process the latter terms, consider the definition of the source ellipticity dispersion (e?e^*) - 5ijcrl. Moreover, as the 
Universe is isotropic, we expect the combination (ej^e^, which has a net orientation, to vanish. From these two complex correlators 
one concludes for the ellipticity components 

K4J = SiAi^Y ' "''^ = fi' 2) ■ (16) 

The source ellipticity terms are processed via this relation and result in o"^/2«, so that they cancel with the last term in (fTSl l. In case 
of E-modes, the remaining shear terms yield 

(PeW) = — ^ r d2^(|ri(f)|2cos22^+lr2(Ol'sin22/? + 2{7i(^)7;(^) + r*(^)f2W}sin2yScos2yS) . (17) 
AAR{t) JAR(t) 

The shears are replaced by means of the expressions 

/^eW = riW cos 2/3 + 72(0 sin 2j8 ; ^b(0 = -ri(0 sin 2y6 + 72(0 cos 2y6 , (18) 

derived from ([8]l. Via (|6]l one then inserts the power spectrum. Formally, at this stage a Dirac delta-distribution with argument 
appears, but as terms quadratic in the ellipticities are involved, the approximation ( fT4] i applies, leading to 

(Pe(^)) = — ^ r dHilnfmPEit)^-^^ r P^{e) ^ P^{t) . (19) 

\ ' AARie) Jam Ar({) Jam 

In complete analogy, one can show to be an unbiased estimator of the power spectrum, too. In Sect. |6] we will drop the assump- 
tions on A(f) made here and examine how well the estimators (fTsT l are suited for more realistic situations. 

4. Covariances 

Next we are going to calculate the covariance of the estimators defined above. Introducing the shorthand notation APx(i) = PxiO - 
Px{i), the covariance reads 

Cov(Px; (, n {APx(£)APx(h} = {Px({)Px(h} - PxC^Pxih ■ (20) 

For simplicity we keep to the case of the E-mode power spectrum in the following; the B-mode calculation is analogous. Plugging 
in the estimators yields for the correlator 

<Pe(W^"')) = r ^TTTTtA -2 ffr^>-. d^' 2/3 + 62(f) sin 2^|' \h (H cos 2/3' + hit) sin l/sf) (21) 
Jam n^AARif) JaM') n^AAR{(') V ' / 

- ^ [ -2'!! r.. (0 cos 2;6 + 62(f) sin 2y3p\ - ^ f -^il^ /|ei(f')cos2;S' + 62(nsin2;S'f\ + ^. 
2n Jam n^AAR(i) V ' I 2n JaM') n^AAR{{') \' / Arf- 
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By comparison with (fTSl l the integrals of the second and third term amount to Pe(£) + cr^/^h and Pe{^') + crl/2h, respectively. 
The expansion of the first term results in a number of four-point correlators of ellipticities. If the shear field is assumed to be 
Gaussian, which is realistic if one considers scales sufficiently large that effects of non-linear density evolution do not have an 
effect on the shear field, this can be written as a sum of products of two-point correlators. The four-point correlators of the intrinsic 
source ellipticities can be decomposed analogously. So terms of type (e„(f)e^(r)) remain, which can be further processed via ( fTSl ) 
as follows: 



2 , 



(22) 



where for the first step we assumed that intrinsic ellipticities and shear are not correlated, while in the second step we used (fTSI l. 
The term in brackets can be recognized as h({ - ('), whereby we can approximate the equation as done before, leading to 

0-2 

{Ump')) = n\%({)y;({')) + 5,p^{2nfh 5^^\{ - t') . (23) 

Correlators of the form {£„{{)£^((')} and {£*(()£^({')} can be dealt with analogously by employing €*({) - £„{-(). 

By expressing the correlators {faif)y^if')} in terms of correlators (X^x(^)'^y(^')) (^^^ X,Y= {E, B]), using (fTsT l, (1211 1 can be written 
in the form 

. - . - r dV r o-l i - - \ o-^ 

with 

Se = ^(2nfn' [_i^(^E(f);f- (^)) + _i^(^E(r)^* (f)) + 4 6^^\{ - r)(^E(0%(0>) 
2 \(27r)^ {2ny j 

= ^(27r)2«3 (-^Pe(^) + T^PEin + 4 6^'\e - t') A Pe(o) , 

where in the second steps, respectively, we have made use of (|6|l together with ( fT4] l. Here and in the following, ^ denotes cosmic 
variance terms, C is the shot noise contribution proportional to cr^ and S stands for the mixed term. Shot noise is caused by the 
intrinsic ellipticity dispersion of the source galaxies and usually dominates on small angular scales (in the Gaussian approximation). 
The fact that a data field of finite extent is used to estimate statistical measures leads to another noise component called cosmic 
variance, which prevails on scales larger than about 5' or { < 2000, respectively, see e.g. iHu & White (.2001,) . 
By subtracting the product of the mean of the estimators, ( l24b turns into 

{APE(e)APE(h}^2-(2nfAn^ [ ,ff - f ff,,J ^\( ' H f« Pe({) + . (25) 

Jam n^AAR({) JA,(h n^AAR{{') \ 2 j 

Note that the Dirac delta-distribution can only have a non-vanishing result if the two integration areas AR(t) and Aj;(^') overlap, i.e. 
if ^ = {'. Performing the integrations and repeating the above considerations in the B-mode case, one arrives at 

/ 2\2 

471 / ,- o-i\ 



{APx({)APxin) = -J- Px{{) + ^ % for X = {E, B] , (26) 
AiAi; \ 2nj 



where the area of an annulus with th i cknes s A(,Ar{£) - 2n(A{, has already been inserted. Equation ( 1261 ) for E-modes is in agreement 
with the results presented in lKaiseii (Il998l) . A similar calculation yields 

<APe(^)APb(^"')) = , (27) 

i.e. as expected the E- and B-mode power spectra are not correlated. 

By means of (|9]l the covariances of all possible combinations of second-order cosmic shear measures can now easily be obtained. 
For instance, one gets for the covariance of the correlation function ^+ 



{AUdx) A^+(02)) ^ ^ Z ^ i'AfumUt'Oi) ((aPe(^) + APM) [aPeH') + AP^it'))) 

1,1' 

" h^Y/ i'^f^^himUl'Bi) ({APE(t}APE{n) + <APb(^)APb(^"'))) , (28) 
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where we have used a discretized version of the transformation because the power spectrum covariance is given for ^-bins. After 
insertion of ( |26] l the covariance depends linearly on A^; therefore, considering the limit Af — » 0, the sum transforms back into an 
integral. 



1 r 

Jo 



Similarly the results for other combinations of second-order measures are derived, so that for the case of correlation functions one 
arrives at 

<A^_(^?i) ^-(m ^ ^^^^ ^ J4(Wi)J4m) i(pE({) + si + IpbH) + ^] 1 , (30) 



2n I \ 2n 

^ ^ _ i 



whereas for the aperture mass we get 

{MM^m) AiM^m)) = At jimvlm) {p^it) + . 02) 

Up to now the covariances are continuous functions of although in reality angular scales are binned. If the binning is sufficiently 
small and the integrands in the above equations are smooth functions of 6i and 6*2, one can simply replace the angles by the central 
values of the bins. However, the shot noise term of ( |29] | and dSOl l is not smooth, but can be further processed by means of the 
orthogonality relation of the Bessel functions 

f de e J^(Wi)J^(W2) = ^ 6(91 - O2) . (33) 
Jo 9i 

The transition to a discrete set of angular bins is done by the relation between Kronecker and Dirac 5's: de^e^AO^^ — > 6^^\02 -61) 
for A6 — > 0. Introducing the number of galaxy pairs in an annulus of bin radius and thickness AO, Np{6) - 2n0A6A n^, one can 
write (|29] l and (|30] | in the directly apphcable form 

<A^i(0i ) A^^(02)} = ^ r die Ay4m)Jo/4m) {Pl({) + Pl({) + ^ (Pe({) + Pb({))] + ^e.e. ■ (34) 

nA Jo [ n ) Np(6i) 

The shot noise term in OTl i cancels, so that in this case no further processing is needed. 



5. Equivalence to real space covariances 

In Paper I covariances of the correlation functions and other second-order measures were obtained by starting directly with an 
estimator for ^± and by taking an ensemble average over all galaxy positions of the resulting covariances. These calculations were 
based on the same assumptions as this work; therefore our results and those from Paper I are equivalent, which is proven in the 
following. 



5.1. Correlation function covariances 

Starting from ( |34] |. the covariance for ^+ can again be split up into cosmic variance, shot noise and mixed term as follows 



1 r°° 

^++ = dee [pl(e) + pl({)}jom)jQ(m 



= dee {PE(e) + PBi{)]Jom)Jom), (35) 

jrAn 



where here and in the following, we drop the arguments di, 02 of the partial correlators etc. The quantity C+ + is already in the 
desired form of Paper I, (27). Inverting the first two equations of (|9]) by means of the orthogonality relation (33[ results in (see also 
S02) 

d6 {^+i0)Jom + ^-{0)JAm] ; /'b(o = ^ d0 {^+{,0)hm - ^-mjAm) , (36) 

Jo 
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from which one concludes Pe(^) + ^'b(^) - J d^9 ^+(9)s'^ ''- If one inserts this expression and the definition of the Bessel function 
in the form Jq(x) - J^'' ^g±u-cosy; ^^^^ fi5[ . the mixed term yields 

^2 n n n^^ r^" 



- 8^4^^ / J d^e j| d^i j| Aip2 e'^<«-«'+«^) ^ j| ( ^02 + g2_20j02Cos^) . (37) 

After integrating over t and the argument of ^+ depends only on the difference ip - (pi - f\ the polar angles ifi, if 2 of 0i, 62, 
respectively, so that one of the i^-integrals can be carried out. The last equality in dJTb is now equivalent to (32) in Paper I. 

Due to the free choice of the sign of the exponential in the Bessel functions Jq and J4 in this and the following calculations the 
results should be invariant with respect to the signs of the angles in the exponent, in this case 6\ and 62- This is indeed the case as 
can be shown by redefinition of the relevant integration variables. We choose the signs of the Bessel function exponentials such that 
the results of Paper I are reproduced in their original form. 

In order to derive the cosmic variance term, we make use of ( [36] l to write 

Plit) + pI{() = ^ J^d^e J" d^6»' (^+(0)^+(6i') + ^_(6»)^_(6»')e'*'^^»"'^'''') e"'-'^"^'^ , (38) 
which then yields, after writing the Bessel functions in exponential form, 



= '^'^ / '^'^ / jl'd^i j|'d^2e'^<«+«'-«'^«^>(f,(0)^,(0')+^-(0)^-(0')e''<^"-^"'>) 

= ^5^JdVJ^ d^i d^2(^+(<Aa)^+(<Ab)+^(<Aa)^-(<Ab)e'"'^-'^^') (39) 

2 r^co f^n r^n r'oo /-»2;r /^27r 

^ kA j d0 0j d(^i J d^2 ^+(iAa)^+(iAb) + 2^ J d(p4>j d(pi J d^2 ^-(tAa)^-('Ab)cos4(^a - (^b) , 



where we have integrated over t and 6' in order to arrive at the second equality and successively defined the vectors <f> :- 6 + 61, 
iff^ := <f~(^i = ^and^b :- ^ + 62 to adjust our notation to that of Paper I. The quantities i^a and i^b denote the polar angles of and 
i//^,, respectively. It is easily verified that (|39] l is identical to (34) of Paper I, so that the equivalence of (A^+(6'i) A^+(6'2)) is proven. 

Incomplete analogy we repeat the calculations for<A^_(0i) A^_(02)) and {A^+{0i) A^_(02))- Note thatC-- =C++ andC+- = 0. 
The mixed terms then read 

= f (ff f d^0 f difii f d^9^+(0)e'''<®-^'+^^^e'*'<^'--^'> = ^ f d,^ ^+ ( Jfl? +0^.- 26162 cos cos 4ip (40) 

SjT'^An J J Jo Jo nAn Jo \ V ' ^ / 



and 



nAh 



- Jd^^ J A^e '^Vi d^2 ^_(0)e'^ («+«'-«^)e4'<^»-^^> (41) 
Aifi^-i^ ^6^ + 6l- 26162 cos ifi^ I i^^^'^y ^1 ^2'^ cos(yV) I (fi 26162 cos ifif^ . 



In the last step we have made use of a geometrical identity based on the cosine theorem, cos(f - (f2) = 
{01 cos(^2 - <Pi) - Oi] I (1^2 - ^il), where ^ denotes the polar angle of 62 - Oi. When replacing P^({) + P^{{) again as in dJST l. 
we find 



-j^ y-»oo /-»2n r*2n 

= ::r~r d(/>(/> d^i d9?2 ^-(iAa)^-(i/'b)cos4(9?i + (^2 - v^a - i^b) 
27rA Jo Jo Jo 

-j^ r^oo /-*27T r*27r 

+ d</></>j^ d^i d^2^+(iAa)^+((Ab)cos4(^i -^2) ■ (42) 
Finally, now writing 

Pli^) - P^i) " 5 /'^^^ / {^AS)^-iO')&''^''^' + ^-{6)^^(6')&^'f) , (43) 



again derived from ( l36l l. we get 

-j^ y^oo /-»27r r^2n 

= — d0 dipi \ d^2 ^+(<Aa)^-('Ab) cos 4(^2 - ^b) . (44) 

Jo Jo Jo 

Hence, the results for all coiTelation function covariances are equivalent to those of the real-space approach of Paper I. 

Moreover, this equivalence can be verified numerically, too. The explicit evaluation of d29] l to ([32]) is straightforward as it merely 
consists of a one-dimensional, though highly oscillatory integral, in contrast to the more complicated formulae derived in Paper I. In 
order to allow for direct comparisons, we make use of the same cosmological model as Paper I, i.e. a ACDM universe with Q.^ - 0.3 




Fig. 1. Covariance of the correlation functions. Left panel: Contour plots for (A^+(6'i) A^+(02)} (thick curves) and (A^_(0i ) A^_(6'2)), 
each with shot noise term removed, which would only yield a bin-size dependent contribution to the diagonal. For {A^+{0i) A^+(02)} 
contours are linearly spaced, starting at 10""^ in the top right corner and with highest value of 10"**. For (A^_(0i) A^_(02)) contours 
are logarithmically spaced, differing by factors of 2.25 each with a maximum value of 10"'*^ in the top right corner. Right panel: 
Contour plot for {A^+(Oi) A^_(02)). Contours are logarithmically spaced, differing by factors of 1 .5 each. Solid contours correspond 
to positive values with maximum at 4.4 x 10"'" in the top left corner, dashed ones to negative values with minimum at -1.5 x 10"'" 
in the top right corner. 



and Qa - 0.7. The power spectrum of density fluctuations is characterized by the primordial slope ns = 1 , the shape parameter 
r = 0.2 1 and the normalization erg = 1 .0. For the computation of the l inear power spe ctrum the fit formula of iBardeen et alj (Il986h 
is used, together with the description of the non-linear evolution by iPeacock & Do dds (1996). To calculate the projected power 
spectrum, a galaxy redshift distribution p{z) oc exp {-(z/zo)'^} with zo - 1.0 and /? = 1.5 is assumed. Finally, the survey properties 

are specified by a galaxy number density « = 30 arcmin"^ and an intrinsic galaxy ellipticity dispersion of cr^ - 0.3. Since all 
covariances scale with the field size as A"', we have chosen a fiducial size of A = 1 deg2 for Figs.[T]and|2] 

In Fig.[T]the covariances of the correlation functions with the shot noise term removed are plotted. We can affirm the general char- 
acteristics of the covariance matrices; however, there are deviations from the results of Paper I exceeding numerical uncertainties. 
For (A^+(0i) Af+(02)) we find an offset of the order 10% which is nearly constant in the plotted region, whereas (A^_(0i) A^_(6'2)) 
quickly tends to zero with increasing separation from the diagonal instead of decreasing further to negative values. Deviations are 
largest for the lower right corner of the (A^+(0i) A^_(02))-plane where we obtain negative values close to the diagonal and values 
close to zero off the diagonal. These differences can be traced back to the fact that in Paper I the correlation functions used as input 
are approximated by power laws for angles larger than 5 deg, which has a non-negligible effect on the results due to the infinite 
range of integration over an gles. Moreover, we find very good agreement between the three covariances of Fig. [T]and those given in 
iKilbinger & Schneideil ('2004). Paper II of this series, who numerically calculated the covariances from the formulae given in Paper 
I without taking the ensemble average over galaxy positions. 

The form of the covariance matrices can be explained in terms of the filter functions introduced in (|9]l. Jo{x) and J4(x) behave 
similarly for large x, but in contrast to J4{x), which tends to for small x, Jo{x) remains positive and is therefore a much broader 
filter. As a consequence, ^+{ff) is strongly correlated with all ^+ at scales 9 < 6', which leads to the nearly square contours in Fig.[T] 
The narrow filter J4{x) results in a quick decorrelation of ^_ for increasingly different angular scales and thus in a diagonal structure 
of the covariance matrix. 



5.2. Aperture mass covariance 

Paper I also gave the covariance of the aperture mass for an estimator Ai{ff) of (M^^), which is a function of the correlation function 
estimators and reads 
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Fig. 2. Covariance of M^p. Left panel: Contour plot for the covariance of M^p. Contours are logarithmically spaced, differing by 
factors of 1 .5 each with highest value of 4.4 x 10 in the lower left corner. Right panel: The thin solid curve depicts the square-root 
of the auto-variance of M^^. In addition, the same measure, calculated from the scaled expression in Schneider et al. (1998), see 
(|48]) . is displayed as thick long-dashed curve. The analytic expression for the minimum auto-variance (dashed line) in absence of 
cosmic variance is plotted, too. {M^ ){9) is shown as thick solid line for comparison. 



where A§ is the angular bin width, the §i are the centers of these bins, and furthermore 6 — mA& with m being an integer. The 
functions T+ and r_ are given in S02. This form of the estimator for (M^p) is insensitive to B-modes. Since the covariance of 
the aperture mass in Paper I is given in terms of the coiTelation function covariances, one can directly plug in the results of the 
previous section and show the equivalence to ( l32b . However, the transformation from ^-t to (M^p) is merely based on (|9]l, so that the 
equivalence proven in the foregoing section readily implies the equivalence also for the aperture mass covariance. 

In absence of gravitational lensing the shot noise term Cap yields the only contribution and is therefore the minimum covariance. 
Either from ( [32] l or from Paper I, (42) one obtains 

0-4 51& 



Note that in the absence of B-modes Cap constitutes the full covariance matrix of {M\). If one considers the auto-variance, the 
integral in ( |46] | can be evaluated analytically giving 



in agreement with the numerical evaluation of Paper I. 

Figure|2]displays the numerical results for the covariance matrix of (M^p) when using the same survey geometry and cosmolog- 

ical parameters as before. We choose again a fiducial field size of A = 1 deg^. Due to the very narrow filter function of (M^p), the 
coiTelation is strongly concentrated towards the diagonal. When the ratio of angular scales exceeds 2, the covariance has decreased 
to less than a third. The structure of the diagonal is analyzed further in the right panel where the square-root of the auto-variance 
is shown. On small scales the noise of intrinsic galaxy ellipticities dominates (at least in the Gaussian approximation), whereas 
on large scales cosmic variance takes over. The overall functional form is identical to the results of Paper I, the differences again 

stemming from their approximation of correlation functions as discussed in the fore going section. 

Moreover, we compare our result for the auto-variance with the one derived in lSchneider et al.l (Il998b who started with a direct 
estimator for (M^p). In the Gaussian case the approximate solution therein reads 



(48) 
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where G is a factor of order unity depending on the choice of the filter function Q for Map; for our choice, G - 6/5. Ni is the 
number of statistically independent apertures which can be placed on the survey area. We set :- rj&^/A, where 77 is a fudge 
factor determined from the arbitrarily chosen condition that the shot-noise terms of both approaches should have the same value. 
The resulting auto-variance is also plotted in Fig. |2] being identical to the power spectrum approach within a few per cent. The 
actual noise of the (M^p)-estimator will be larger though, because 77 ~ 0.6 is smaller than unity. This is expected because for the 
estimators ( fT5] l every ellipticity in the field is used, whereas in placing circular apertures information is left out even if one allows 
the apertures to partly overlap. 



6. Generalizations 

6.1. Shear tomography 

The inclusion of photometric redshift information about the distance of the observed galaxies allows one to measure the effects of 
gravitational lensing of the large-sc ale s t ructur e in three dimension s, enabling the determination of cosmological parameters with 
considerably tighter constraints, e.g. i Hul (Il999l) : iBacon et"aLl (l2005 l). In particul ar, it is pos sible to follow the evolu tion of structure 
with time, which depends on the density of dark energy and its equation of state dHu & Jain 2004; Benabe d & van Waerbeke 2004- 
Besides, redshift inf ormation is essential to eliminate s ystemat ic effects such as intrinsic alignment of galaxies from the cosmic 
shear signal, see e.g. iKing & Schneider (.2003); Hirata & Seliakl (|2004.). Therefore it is of interest to review the results obtained for 
the power spectrum covariance in the context of shear tomography. Note that we retain our assumptions used in Sects. [3]to|5] 

Consider A^bin redshift bins into which all galaxies can be sorted. For each bin k an equivalent lens plane with convergence /Cg ' 
as well as ellipticities e*^*^ and shears y**' are defined. Then one can construct auto- and cross-correlation tomography power spectra 

= (27Tf6^^\e - n Ps\{) , (49) 

{4Wt(n) = , 

where we again assumed the E- and B-mode cross power spectra to vanish due to parity symmetry, i.e. P^^^it) = 0. The tomography 
power spectra are connected to the three-dimensional matter power spectrum via ^ with modified lensing efficiencies g'^^'Of) = 
j^'' dx' p'''^''(x')(^ -xlx') for bin k. Here p'^^^ix) stands for the distribution of galaxies within bin k. The distributions of different bins 
will in general overlap, for instance due to uncertainties in the determination of photometric redshifts. From the power spectra one 
obtains second-order shear tomography measures and their interrelations, as for instance 

^f\e) = (rfr®) (0) ± (r«r?) (») - [ |^ J^m [pfn) ± pfio] . (50) 

Let the number of galaxy pairs in bin k be A^**^\ and the corresponding number density n^*^' - N^'^^/A. Then (l23T l is generaUzed to 

(e«(f)f *(r)) = {t\Off\r)) + c5„,5,,^(2;r)2^W 6^^\{ - {') . (51) 

The 6ki in the second term appears because intrinsic ellipticities of galaxies are only correlated with themselves, so that there cannot 
be a correlation between different bins. Now we define an estimator for the shear tomography power spectrum 

Pfii) ^(k)-jAAR{() L (^f'W'^°^2yS + e^*'(f)sin2^)(e®*(^)cos2/J-H4'>*(f)sin2^)-5„^ , (52) 

and likewise for P^^'\i). Both estimators are unbiased, i.e. {P^"\i)) = P^^Hl) for X = [E, B] and all k, I. Repeating the calculation 
as in Sect. 21 one arrives at the covariances of the power spectra, which read 

{APfmpfii')) = [pfmfCi) + Pf(()Pf\{)) %, with -Pf := Pf + . (53) 

This r esult is in agreement with the covariance derived for instance in iTakada & JainI (l2004l) . iMa et alJ (l2006l) and iHuterer et al.l 
(1200 6'). Note that as in the non-tomographic case, any E- and B-mode cross-covariances of the tomography power spectra vanish. 
By means of the second equality in (ISOl l and corresponding relations, the covariances of shear tomography second-order measures 
can now easily be obtained. We confine ourselves to the example of for which one gets in the case of P^(() = 

1 r°° r _ _ _ _ cr^ 1 

{A^:'\ei) A^f(02)) ^2^j^ ^ Jom)Jom) Upf\{)P^f({) + Pf({)P^f(i)) + -^r^r^ {dikdj, + 6u6jM . (54) 

ISimon et al.l ([2004") calculated such a covariance matrix from mock galaxy catalogs drawn from Gaussian random fields by means 
of a Monte Carlo simulation. As a practical ex ample we rec alculate this matrix, the result given in Fig. [3] We employ the same 
cosmological model and redshift distribution as ISimon et al.l (2004), which are identical to the ones used in Sect.|5] Similarly we 
adopt a mean galaxy density of n = 30 arcmin"^, an intrinsic ellipticity dispersion of cr^ - 0.3, a single square field of observation 
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Fig. 3. Covariance matrix 
of tomographic shear 
correlation functions 
plotted without shot 
noise term. Each block 
consists of the covariance 
matrix of the correlation 
functions given on the 
margin. Within a block, 
entries correspond to 65 
logarithmically spaced 6- 
bins in the range between 
2' and 40'. Contours are 
logarithmically spaced 
with the lowest values of 
about 10"'° in the upper 
left corner (darkest shad- 
ing) and largest values 
of 2 X 10"^ in the lower 
right corner (brightest 
shading). 




of size 1 .25 X 1 .25 deg^, and four equally sized disjunct redshift bins between and 3. Angular scales range from 2' to 40', divided 
into 65 bins. Again we leave out the shot noise term in the plot as it would only contribute to the diagonal of each sub-matrix. 

The overall behavior of values in Fig.[3]is the same as in lSimon et al.l (|2004|) . the apparent differences in the structure within the 
blocks being due to the linear spacing of the 0-bins in lSimon et al.l (|2004|) . As expected, the structure of blocks of type ^A^+'^ '^^+'^) 
is similar to Fig. [T] The large values for the covariance of the correlation functions at the highest redshifts are caused by both 
increasing P^^\() and the low number of galaxies in the highest redshift bin. 



6.2. Arbitrary survey geometry 

Data of real cosmic shear surveys will in general not satisfy our assumption of a simple, filled geometry, but contain holes and gaps 
due to the masking of stars, pixel damages or gaps between CCDs. Moreover, it has to be checked whether the assumption that all 
relevant angles are much smaller than the survey size is well justified for realistic configurations. In order to attain results of a more 
practical value, we now take no priors on the survey geometry, but retain the assumption of randomly distributed galaxy positions, 
and repeat the calculations of Sects.[3]and|4] 

A{{) is now used in its general form; it can be calculated at least numerically for realistic cases. The estimators ( fT5] l for the power 
spectra are kept, but will not be unbiased any more, so that we first recalculate their expectation values. Splitting the ellipticities 
into shear and source ellipticity terms via ( fTTT l. one arrives in case of the E-mode estimator at 



AARie) Ja, j j 



(55) 



: (ri(^i)rI(^2)cos2 2/3 + r2(/'i)r;(^2)sin2 2/3 + yi{ti)y2m sin2yScos2yS + nihWiih) sin2yScos2/?) 
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While the terms containing source ellipticities are processed as before and cancel with the cr^/2n term, the shear coiTelators can be 
written in terms of convergences by means of the real and imaginary part of ([8]l as for example 



(fi imiifi)) = (MhKifi)) COS 2/3i cos 2/32 + {MhTKlih)) sin sin 2^2 
- {Mh)T<lih)) cos ipi sin 2^2 - {MhrKlih)) sin ipi cos 2^82 



(56) 



where denotes the polar angle of f,-. One proceeds likewise for the other shear coiTelators. After replacing the coiTelators of the 
convergences by power spectra via (|6]l, one obtains an expectation value 



Ar{() Ja, 



JAnit) 

where we have defined 



= ^ J d^A A(/' - A)A*(r - /"i) {^'e(^i) cos 20Si-^) cos 208i-/?') + ^b(A) sin 20Si-yS) sin 2(/?i -/?')) - 
Repeating this calculation for results in an analog to ( fSTl l. now with 

<Un = ^ J A^€y A(/'-/'i)A*(r-A){^B(A)cos2(/?i-yS)cos20Si-/3') + PE(A)sin20Si-yS)sin206i-^')} ■ 



(57) 



(58) 



(59) 



In the special case that the assumptions of the previous sections hold, the observed power spectra are equal to the theoretical power 
spectrum. However, in the general case they are a convolution of P^{{) and P^{{) with the power spectr um of the field aperture , 
modified with factors stemming from the E/B-mode decomposition. This is a well-known result, derived by'Baumg art & Fry! (Il99lh 
in the context of galaxy redshift surveys for arbitrary weighting functions. Similarly, we can allow 11(0) not to just have values of 
and 1, but assign to every point of the survey an individual weighting. The weighting fu nction, i.e. in our cas e HiO), can then be 
optimized such that the variance becomes minimal. For the case of galaxy redshift surveys, jpeldman et al.l(ll994l) derived a minima l 
variance weighting under the assumption of Gaussian fluctuations; a more gener al consideration can be found in lHamiltonl(ll997h . 
This optimization could be used as an alternative to the numerical approach by iKilbinger & Schneided (|2004|) . in order to design 
optimal survey geometries. 

The covariances of the estimators are calculated in the same manner as in Sect.|4] After some algebra one obtains for X - [E, B] 



{APx({)APx(h) = ^ ^-^ f r AH' 2 

\ ' Ar{{)Ar{(') Jam JA.Cn 



In 



(60) 



where cr'^{(, (') = [crlilTiflA] A({ - {') cos 2(j3-/3') has been defined for convenience. Now the covariance matrices are in general 
no longer diagonal. Furthermore, a correlation between Pe and Pb is introduced: 



{APE(e)APB((')) = ^ — - r r d^f 2 

Ar({)Ar({') Jam ■JAM) 



2h 



(61) 



where we defined 



Pf;{tJ') = ^ J AHyA({ - h)iA*{t' - h){Pv.{ey)co^2iJ3y -/3)sm2(J3i - /3') - Psifi) sm2(fii -/3)cos2(J3i -yS')} (62) 

and d-l^i{, (') = [(rli2jT)^ /a] A{{ - {') sin 2(J5 - /3'). Note that in the limit A(0 ^ ^^^^(0, the fore going results for a simple survey 
geometry are regained also for the covariances. 



6.3. Some specific examples 

In order to demonstrate the influence of field geometry on our estimators, we consider a number of examples for which A{{) can 
be calculated analytically. First, we assume simple square fields of size 1x1 deg^ and 5x5 deg^ each, as well as a lOOOdeg^ 
square geometry, its size being of the order of currently planned cosmic shear surveys. Second, we place a hole of 0.25 x 0.25 deg^ 
at the center of the 1 deg survey. Finally, a chessboard configuration is examined, where 8 sub-fields of size 0.125 deg^ each are 
embedded in a 2 deg^ field to get the same effective area as for the 1 deg^ survey in the first example. By letting A^ — » in (l57T i. 

we calculate a phase-averaged observed power spectrum, {P^'^)i{) = 1 /27t A/3 P J*(f , {), for each geometry. We set Psii) = 0; 

however, does not vanish since it has contributions from Pe{{), see ( fSST l. 

In Fig.|4l left upper panel, the resulting phase-averaged observed E-mode power spectra for all field configurations are plotted. 
In addition, the theoretical power spectrum as well as the observed B-mode power spectrum in the case of the 1 deg^ square field are 
depicted. The lower left panel shows the amplitude of the observed E-mode power spectra relative to the theoretical one. Moreover, 
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Fig. 4. Influence of survey geometry on the estimated 
power spectrum. Upper left panel: E-mode power spec- 
tra {P^'^){{), averaged over the polar angle of t, ob- 
tained for different survey geometries. The sparsely 
dotted curve shows the original E-mode power spec- 
trum PE(f)- The solid line corresponds to a square field 
with 1x1 deg^, the dotted one to a square field with 
5x5 deg^ and the long-dashed one to a square survey 
with 1000 deg^. In the case of the short-dashed curve 
the 1 deg^-field contains a hole, while the dot-long- 
dashed line corresponds to the chessboard configura- 
tion. For comparison, the corresponding observed B- 
mode power spectrum (caused by E-modes) for the 
1x1 deg^ square field geometry is plotted as dot- short- 
dashed curve. Lower left panel: Fraction of observed 
E-mode power spectrum over original E-mode power 
spectrum. The coding is identical to the left upper 
panel. Right panel: Power spectrum of the phase- 
averaged field aperture {A^}(£). The coding of the 
curves is the same as in the left panels, the curve for the 
chessboard configuration lying on top of the solid line 
in the central part of the diagram. The axes are normal- 
ized such that <A^)(0) = 1 for the 1 deg^-field. Note that 
the maximum of the dotted curve is at 25. The even nar- 
rower curve for the 1000 deg^ field has not been plotted 
for reasons of clarity. 



the aperture power spectra {A^}({), again averaged over /3, for the different geometries are given in the right panel, for our examples 
all being composed of sinc^-functions. Note that our assumptions about the survey geometry made in the previous sections would 
lead to a Dirac-delta spike at ^ = in this plot. 

For large { all observed power spectra trace the theoretical one well, whereas the configurations with small survey area deviate 
significantly for £ < 1000. Due to the limited field extent no ellipticity pairs are measured at large separations on the sky, corre- 
sponding to a lack of information about the power at small £. In the case of larger effective survey size, deviations are therefore 
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smaller and start at smaller angular frequencies. The relative amplitude of the observed power spectrum decreases below 90 % at 
t < 300 for a survey size of 5 X 5 deg^ and at ^ < 50 for the 1000 deg^ survey. The convolution kernel of the 5x5 deg^-field is 
already quite narrow, the Dirac-delta distribution therefore being a reasonable approximation in the case of large survey areas. 

Cutting a hole into the field has very little effect on the observed power as there is still information available on all scales; only 
the noise is expected to increase. The relative height of the side lobes of the aperture power spectrum increases moderately, so that 
{P'^^){t) is slightly lowered. On the contrary, the curve for the chessboard configuration lies significantly lower than for the 1 deg^ 
field. This can be explained by the pronounced side lobes of the otherwise identical aperture power spectrum, which are generated 
by the periodicity of the field geometry. Power is already reduced considerably on angular scales corresponding to the size of the 
chessboard squares, leading to an early deviation from the theoretical power spectrum. 

As expected, the observed B-mode power spectrum is significantly smaller than the E-mode signals for larger angular frequen- 
cies, but for { < 200, corresponding to angular scales which do not fit into the field, the B-modes attain the same order of magnitude 
as the observed E-modes because in this range the phase-averaging of the sin^ 2(fii - fS) and the cos^ 2(fii - p) term in yields 
the same result or, illustratively, E- and B-modes cannot be distinguished anymore on large angular scales. This means that apart 
from a bias one also has to take into account a leakage between E- and B-modes, which is severe in case one would like to use the 
power spectrum estimator ( fTsT i to measure cosmological B-modes. Since these are potentially significantly smaller than the E-mode 
signal, even a small leakage will then cause a strong contamination. 

Taking into account the results for these exemplary field geometries, it is unfortunately not so obvious how the calculation of 
the covariances of other second-order measures should proceed in the general case, using dMT l. since ^ has to hold also for the 
estimators, which is not the case if Pe{?) is biased. However, the survey geometries considered here are quite artificial; in the real 
world, one will probably be faced with a large contiguous field, with certain regions, located quasi randomly, been cut out. In this 
situation, the covariance is expected to behave like in the connected field geometry case, with the survey area A taken to be the 
effective area, i.e. with the cut-out regions subtracted - our case with the hole supports this conjecture. The results for the specific 
examples discussed above may then be used to assess the accuracy of this approximation. 

7. Conclusions 

In this paper we presented a new approach to calculate analytically covariance matrices of second-order cosmic shear measures 
under the assumption of Gaussian density fluctuations, uniformly distributed galaxies and a simple survey geometry. By introducing 
estimators for the E- and B-mode power spectra and performing the derivation in Fourier space, we arrived at compact results which 
will simplify further analytical considerations and speed up numerical work. These results are fully equivalent to those found in a 
real-space approach by Paper I. 

We then formulate d the covariance matrices in the context of shear tomography. Compared to Monte Carlo simulations employed 
bv lSimon et al.l ( |2004|) . the analytical ansatz leads to an enormous gain in computational speed. However, the former method has 
the advantage that it is directly applicable to arbitrary survey geometries. 

We also generalized our findings to arbitrary survey geometries and examined the bias thereby introduced on our estimators. In 
this context one might pose the question why the convergence power spectrum, which contains all second-order information about 
cosmic shear, is not directly estimated from the data and solely used for further analysis. However, the power spectrum estimators 
incorporate a Fourier transformation for which galaxy ellipticities have to be measured in principle on the whole sky. This problem 
is avoided by the assumption \ff^ <s: A, but will lead to a bias in case of finite and more complex field geometries. In principle, the 
power spectrum could also be calculated from real-space measures via the inverse of (|9]l, but this is hindered by the infinite range 
of integration over the angular separation 0. 

Although the theoretical power spectrum cannot be inferred to satisfactory accuracy from the estimators (fTSl ). they may still be 
of use in cosmic shear data analysis. As demonstrated in Sect. |6l it is relatively straightforward to calculate expectation values of 
observed power spectra P"*"^ from the theoretical Px( and success i vely c ompare these to results from observations. Another way to 
directly extract the power spectrum was discussed in lHu & White! (1200 lb who determine band powers with a maximum likelihood 
method, which is able to recover the power spectrum irrespective of the survey geometry, but is not suited for large fields due to its 
computational requirements. Furthermore, due to its narrow filter, M^^ can yield a fairly accurate estimate of the power spectrum 
localized at scales { ^ 4.25/0 (Bartelman n & Schneideij|l999 ). In order to determine cosmological parameters though, one would 
not make use of the power spectrum. For this purpose the correlation functions are the second-order measures of choice since they 
remain unbiased for arbitrary survey geometries. 

All methods discussed here are limited by the assumption of Gaussianity, that will break down on small, non-linear scales. 
Anyway, in case of non-Gaussian density fluctuations the power spect rum is not sufficient for their c har acterization any more, s o 
higher-order statistics have to be taken into account. In fact, as shown bv lKilbinger & Schneided (l2005h and lSemboloni et aTl (l2007h . 
the non-Gaussianity modifies the covariance substantially on angular scales below ~ 15', depending on the source redshift distri- 
bution. The increase of the covariance due to non-Gaussianity is so str ong that the shape noise term is nearly always sub-dominant 
compared to the cosmic variance contribution. Sem boloni et aTl (12007^ have used numerical simulations to quantify the deviation of 
the covariance matrix of the shear correlation function from the Gaussian case, giving analytic fit functions. We hypothesize that a 
similar approach is possible in the context of our formulation, where now the fitted correction function for non-Gaussianity will be 
be applied to the power spectrum covariance ( |26l ). We will consider this in future work. 

For the full treatment of four-point correlators needed for the covariances of second-order statistics in the non-Gaussian regime, 
only few analytical methods exist, so that one has to rely on numerical work, for example ray-tracing simulations. As also the 
analytical fit functions for the input power spectra (Peacock & Dodds 1996; Smith et al. 2003) have limited accuracy, simulations 
will be indispensable for the determination of covariances for future cosmic shear surveys, designed to serve as high-precision tools 
for cosmology. 
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Moreover, the calculation of covariances via a power spectrum approach has revealed an interesting analogy to the analysis of the 
Cosmic Microwave Background. Similarly, power spectrum estimators for the E- and B-mode polarization of the CMB are defined 
in an idealized survey geometry, namely on the full sky. Due to incomplete sky coverage or masking of the Milky Way in real data 
these estimators become biased and mix E- and B-modes. However, in t he con text of CMB pola r ization there exist me t hods which 
allow to debias estimators and single out B-modes, see e.g. lHivon et al. I (|2002|) and lSmithl(l2006l) : ISmith & Zaldarriagal(l2006j) . The 
further investigation of these methods, transfered to cosmic shear, may prove valuable. 
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